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Reference:  MAC/ 1 73/01 3 
1.  SUMMARY 


The  visible  radiation  from  a  thermally  ionized  high  density  plasma 
produced  by  explosively  generated  shock  waves  has  been  investigated  as  a 
possible  laser  pumping  source.  It  is  considered  that  the  major  problem  is  the 
efficiency  of  conversion  from  chemical  energy  to  radiation  in  the  absorption 
band  of  the  laser.  The  study  has  therefore  been  concentrated  on  the  mechanism 
of  the  emissive  process  with  a  view  to  increasing  the  radiation  emitted  in  a 
given  spectral  region.  The  real  gas  thermodynamic  variables  behind  the  shock 
waves  have  been  evaluated  for  argon  and  helium  and  used  tc  calculate  linear 
absorption  coefficients  based  on  Unsflld-Kramers '  theory.  Using  Kirchhoff's 
law,  the  radiant  intensities  at  the  surfaces  of  optically  thick  argon  and 
helium  plasmas  at  20,0C0°K  have  been  calculated  and  compared  with  observed 
values  in  the  visible  spectral  region.  The  effects  of  various  parameters  such 
as  gas  composition,  pressure  and  shock  velocity  on  the  spectral  radiance  have 
been  studied  experimentally.  It  is  concluded  that  the  radiation,  from 
explosively  shocked  noble  gases  is  essentially  black  body  in  nature. 


2.  INTRODUCTION 


In  recent  years  there  has  been  considerable  interest  in  the  use  of 
lasers  for  producing  transient  high  radiation  densities.  At  present,  high 
energy  Imer  systems  are  relatively  inefficient  and  very  large  input  energies 
are  needed  if  t  gh  output  energies  are  required.  Also,  since  the  available 
pumping  times  for  present  day  laser  materials  are  short,  this  energy  must  be 
supplied  at  a  very  high  rate. 

High  explosives  are  probably  the  most  compact,  cheapest  and  lightest 
energy  sources  available;  the  explosive  60/40  RDX/ENT  contains  about 
0.77  x  10s  joules  kg-1*  To  store  the  equivalent  energy,  a  capacitor  bank 
would  have  to  be  at  least  50,000  times  as  large  in  both  volume  and  weight. 
However,  since  explosives  are  destructive  in  nature  they  are  useful  only  in 
certain  classes  of  application. 

The  use  of  explosive  energy  to  pump  lasers  poses  the  problem  of  converting 
the  chemical  energy  of  the  explosive  into  radiation  of  the  appropriate 
frequency.  Of  the  several  methods  available  only  one  is  considered  here;  an 
explosively  generated  shock  wave  in  a  monatomic  gas  is  intensely  luminous  and 
the  light  may  be  used  directly  to  pump  lasers.  This  report  describes  a  3tudy 
of  the  mechanism  by  which  the  energy  of  the  expanding  detonation  products  is 
converted  into  radiation. 
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As  a  starting  point  in  this  investigation,  the  argon  flash  bomb  was 
evaluated*  This  is  a  simple  device  designed  to  meet  the  needs  of  high  speed 
photography  and  consists  of  a  cardboard  box  containing  a  thin  slab  of 
explosive  at  one  end  and  fitted  with  a  cellophane  window  at  the  other  end. 
Y/hen  the  air  is  flushed  out  with  argon  and  the  explosive  is  detonated  several 
hundred  million  candlepower  is  produced  for  about  1CT4  see.  Assuming  for  the 
moment  that  the  radiative  intensity  corresponds  to  that  of  a  black  body  at 
20,C00°K  the  total  energy  radiated  is  about  9  x  10s  watts  cm**2  for  1CT4  sec* 
or  90  joules  cm”2  of  radiating  surface.  Considering  a  flash  boiab  containing 
a  slab  of  RDX/TNT  1  cm  thick  the  chemical  energy  liberated  is  about  6.4  *  103 
joules  cm*2  of  explosive.  This  leads  to  an  overall  conversion  efficiency  of 
i«4  per  cent  from  stored  chemical  energy  to  total  radiated  energy. 


3.-  THEORETICAL  CONSIDERATIONS 


In  this  Section,  the  theoretical  aspects  of  the  energy  conversion  process 
are  studied  with  a  view  to  predicting  the  optimum  conditions  necessary  for 
obtaining  maximum  radiance,  in  a  given  frequency  interval,  for  a  given 
explosive/gas  system.  The  emphasis  is  placed  more  on  obtaining  maximum 
spectral  intensity  than  on  conversion  efficiency  from  chemical  to  radiant 
energy. 


3.1  Equilibrium  Shock  Y7ave  Calculations 


An  essential  prerequisite  to  the  understanding  of  the  radiative  properties 
of  a  shock-heated  gas  is  a  knowledge  of  the  equilibrium  real  gas  thermodynamic 
variables  behind  the  shock  front.  Under  the  extreme  conditions  produced  by 
high  explosives,  deviations  from  perfect  gas  behaviour  may  result  from  thermal 
excitation  and  ionization  of  the  gas  behind  the  shock  front.  These  difficulties 
can  be  overcome  by  the  introduction  of  quantum  statistical  mechanics  of  the 
microcanonical  assembly  into  the  shock  wave  calculations  to  allow  for  these 
real  gas  effects. 

There  have  been  many  such  calculations  for  plasmas  produced  by  shock 
waves,  particularly  for  low  pressures  (l),  but  little  attention  has  been  given 
to  conditions  of  very  high  pressure  and  temperature  behind  the  shock  front. 

Bond  (?)  has  calculated  the  equilibrium  properties  of  argon  for  ambient 
pressures  of  10  torr  and  593  torr  for  incident  shock  velocities  of  between 
3  x  105  cm  sec-1  and  9  x  10s  cm  sec”1.  Lowering  of  the  first  ionization 
potential  was  not  taken  into  account  in  his  calculations  and  at  the  higher 
pressures  and  velocities  considered  by  him  this  could  lead  to  considerable 
error.  Second  ionization  and  excitation  were  also  neglected.  Since  accurate 
values  of  temperature,  degree  of  ionization  and  particularly  electron  density 
are  essential  if  meaningful  absorption  coefficients  are  to  be  evaluated,  the 
physical  bases  of  existing  equilibrium  shock  wave  calculations  were  reviewed. 


/Considerable 
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Considerable  progress  has  been  achieved  in  this  field  recently,  particularly 
with  regard  to  the  developeent  of  a  consistent  physical  code!  (l ,3)*  The 
results  of  preliminary  calculations  for  argon  and  helium  are  reported  here. 

It  should  be  enphasised  that  these  calculations  assure  that  equilibrium 
conditions  maintain  behind  the  shock  front.  The  problem  of  energy  relaxation 
is  not  considered,  the  assumption  being  cade  that  the  equilibrium  state  is  net 
influenced  by  the  physical  processes  that  govern  the  approach  to  equilibriue. 


For  a  conatceic  gas,  the  equilibrium  conditions  are  determined  by  the 
equations  for  the  conservation  of  sass-flon,  momentum  ami  total  energy  across 
the  shock  together  with  the  equation  of  state;  the  Saha  equation,  and  the 
equation  for  the  ionization  potential. 


Using  the  flour  regimes  shotm  in  Fig.  i ,  based  on  a  stationary  shock,  and 
assuming  no  ionization  in  the  undisturbed  gas  ahead  of  the  shock  front,  the 
conservation  equations  can  be  written  in  the  fora* : 


Vi 
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where  the  subscript  o  refers  to  conditions  in  the  ^disturbed  ga3  and  the 
subscript  1  to  the  equilibrium  conditions  behind  the  shock  front. 


The  equations  of  state  ere 


D 
'  O 


p  BT 
ro  o 
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and 


?i  =  (l  +  at)  Pt  2Ti 


5 


where  R  is  tno  gas  constant  per  unit  mass.  The  Saha  equation  can  be  written 


/in 


To  achieve  some  economy  in  presentation  the  symbols  in  ccssncn  use  are  not 
defined  in  the  text  but  all  definitions  are  tabulated  in  the  Appendix. 
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in  the  form: 


The  ionization  potential  x  is  equal  to  the  ionization  potential  xo  for  the 
isolated  atom,  reduced  by  AX  due  to  the  surrounding  plasma.  The  equation  for 
the  lowering  of  the  ionization  potential  used  in  the  preliminary  calculations 
was  that  favoured  by  Pomerantz  (4),  who  critically  reviewed  the  literature  on 
this  topic.  This  equation,  which  was  derived  by  Ecker  and  Y/eizel  (5)>  is 
given  in  the  form: 

1  /N 

AX  =  0.67  x  10"6  Nj  +  0.37  x  10“7  (y 


The  first  term  on  the  r.h.s.  which  is  generally  the  dominant  term  represents 
the  ionization  potential  lowering  due  to  the  plasma  microfield  (the  Unsold  term) , 
whilst  the  second  terra  accounts  for  the  energy  level  lowering  effects  arising 
from  plasma  polarization.  This  theory  has  recently  been  criticized  on  the 
grounds  that  it  leads  to  too  large  a  value  of  A£  for  the  densities  involved  in 
most  laboratory  plasmas  (6,7) • 

The  electron  density  is  given  by: 


ci  Pi 

N  =  ■  — 

e  "  1  -  01  KP 


8 


The  calculation  of  electronic  partition  functions  is  an  astrophysical 
problem  of  long  standing.  The  partition  function  is  defined  by  the 
relationship 

Q  =  nlo  gn  exp(”kf)  . - 

where  gn  is  the  statistical  weight  or  degeneracy  of  an  atomic  energy  level  whose 
energy  is  e^.  The  summation  is  formally  over  all  available  energy  levels. 

There  are  two  main  problems  involved  in  the  calculation  of  the  partition  function. 
The  first  is  the  fundamental  problem  of  its  divergence.  This  arises  because 
there  jxists  for  a  given  atom  an  infinite  number  of  discrete  electronic  states 
converging  to  a  continuum,  and  since  the  excitation  energies  of  these  states 
are  finite  the  partition  function  diverges.  This  problem  is  generally  overcome 
by  employing  some  kind  of  cut-off  procedure  whereby  the  summation  of  the  energy 
levels  is  terminated  at  a  specified  value  of  excitation  energy.  In  some  cases, 

/the  ..... 
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the  partition  functions  are  terminated  at  a  value  determined  by  the  equation 
for  the  lowering  of  the  ionization  potential.  The  second  problem,  which 
arises  from  the  lack  of  available  information  on  the  electronic  energy  levels 
of  the  atoms  and  ions,  is  usually  ignored  and  only  the  available  tabulated 
energy  levels  are  included  in  the  summation  (8).  Recently,  however,  McChesney 
(9)  has  estimated  the  missing  energy  levels  and  evaluated  their  effect  on  the 
calculated  equilibrium  thermodynamic  properties  of  shock-heated  xenon. 

3»2  Shock  Nave  Calculations  for  Argon  and  Helium 


The  equilibrium  properties  of  explosively  shocked  argon  and  helium, 
initially  at  atmospheric  pressure,  have  been  calculated.  In  these  calculations, 
no  attempt  was  made  to  compute  exact  partition  functions.  In  the  case  of 
argon,  the  generally  used  approximation  Qq  =  1  and  Q®  =6  was  used,  while 
with  helium,  due  to  the  very  high  excitation  energies  of  the  first  excited 
states  (l 9*8  eV),  the  partition  functions  were  approximated  to  the  statistical 
weights  of  the  ground  states  of  the  atoms  and  first  ions. 

In  these  calculations,  the  following  assumptions  were  made: 

(i)  the  flow  phenomena  are  adequately  described  by  hydrodynamic  theory 
for  a  one-dimensional  shock  wave  with  heat  conduction,  radiative 
cooling,  and  viscosity  neglected; 

(ii)  thermal  equilibrium  exists  in  finite  regions  behind  the  shock 
front; 

(iii)  the  hot  gas  behind  the  shock  front  consists  of  a  perfect  gas  mixture 
of  atoms,  electrons,  and  first  ions;  and 

(iv)  reflection  of  the  shock  wave  occurs  at  a  perfectly  rigid  wall. 

In  the  preliminary  calculations,  multiple  ionization  and  electronic 
excitation  of  the  atoms  and  ions  were  ignored.  This  should  not  lead  to 
appreciable  error  in  the  case  of  helium.  However,  due  to  the  lower  excitation 
and  ionization  energies  of  the  argon  atom,  there  will  be  appreciable  second 
ionization  and  electronic  excitation  which  could  lead  to  error,  particularly 
in  the  calculated  temperature  and  electron  density. 

Since  the  conditions  in  the  undisturbed  gas  ahead  of  the  shook  are  known, 
and  the  shock  velocity  uQ  is  easily  measured.  Equations  1_  to  6  can  be  solved 
by  an  iterative  procedure.  Initial  estimates  for  the  pressure  and  degree  of 
ionization  behind  the  shock  wave  were  based  on  the  assumptions  that  half  the 
energy  behind  the  shock  front  goes  into  ionization  in  the  case  of  argon  and 
that  the  majority  of  the  energy  in  front  of  the  stationary  shock  is  kinetic 
energy.  This  leads  to  the  approximate  forms  of  Equations  2  and  3,  viz. 
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P  =  pu2  . 10 

1  rO  0  — 


and  ai  tf-  =  -r  .....  11, 

MA  4 

With  these  initial  estimates  Equations  1.  to  6_  were  solved  with  x  =  X0 • 

From  the  approximate  values  of  ai ,  Pi  and  Ti  thus  calculated.  Equations  7  and 
8  are  used  to  determine  the  first  approximation  of  x»  Equations  1_  to  6  are 
then  resolved  with  the  new  value  of  x •  This  process  is  repeated  until  the 
variables  converge  to  the  desired  degree  of  accuracy. 

The  shock  equations  were  solved  for  a  shock  wave  propagating  at 
8.3  x  1^5  cm  sec"1  (Mach  27)  into  atmospheric  pressure  argon.  After  about 
twenty  iterations,  the  following  values  were  obtained,  Ti  =  19,880°K, 
ai  =  0.75,  Pi  -  1*39  x  10"z  g  cm"3.  Pi  =  1.009  x  109  dyn  cm"2  and  N  = 

1.57  x  102°  cm"3.  e 

Equations  1_  to  were  also  solved  for  atmospheric  pressure  helium  with  an 
incident  shock  velocity  of  10.1  x  105  cm  sec"1.  Since  the  Mach  number 
(about  10)  and  degree  of  ionization  are  low,  no  lowering  of  the  ionization 
potential  was  considered.  In  fact,  the  calculated  thermodynamic  variables 
differed  only  slightly  from  those  obtained  using  the  usual  Rankine-Hugoniot 
equations  for  an  ideal  gas.  The  calculated  values  were  Ti  =  9465°K, 
ai  =  2.7  x  10-6,  Pi  =  1.66  x  10"4  g  cm"3,  Pi  =  1.27  x  10s  dyn  cm-2  and  N  = 

2.65  x  101  6  cm"3. 

In  the  case  of  reflected  shocks,  an  additional  equation  is  needed  to 
describe  the  boundary  conditions  imposed  by  the  reflecting  surface,  i.e.  the 
gas  particles  behind  the  reflected  shock  are  rendered  stationary.  In  a 
coordinate  system  moving  with  the  reflected  shock,  this  equation  can  be 
written 


*oR 


-  u 


1R 


=  u 


-  u. 


12 


where  uq^  and  u^  are  the  flow  velocities  of  the  gas  entering  and  leaving  the 
reflected  shock. 

With  the  aid  of  Equations  1_  to  6  and  the  boundary  condition  described  by 
Equation  12,  the  physical  variables  behind  the  reflected  shock  in  helium  were 
obtained.  They  were  T  =  20,370°R,  a  =  9*6  x  10"3,  p  =  6.45  x  10~4  g  cm"3, 

P  =  1.009  x  109  dyn  cm"2  and  N  =  2.43  x  101  8  cm.  No  consideration  has  yet 
been  given  to  the  problem  of  the  lowering  of  the  ionization  potential.  The 
inclusion  of  ionization  potential  lowering  effects  would  result  in  a  somewhat 
higher  degree  of  ionization  and  electron  density  and  a  lower  temperature. 

/Equilibrium  ..... 
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Equilibrium  shock  wave  calculations  are  in  progress  which  allow  fcr 
multiple  ionization  and  excitation  and  the  lowering  of  the  ionization 
potential.  These  effects  are  likely  to  be  of  considerable  importance  in,  for 
example,  the  reflected  shock  in  argon.  The  Debye-Huckel  theory  of  point 
charges  as  applied  by  Griem  (?)  for  multiple  ionized  plasmas,  is  being  used  as 
the  basis  of  these  calculations.  This  theory  provides  an  internally  consistent 
physical  model  for  lowering  of  the  ionization  potential  and  terminating  the 
partition  function. 

3«3  Radiant  Emission  from  the  Hot  Dense  Plasma 

The  energy  radiated  from  a  hot  dense  plasma  arises  from  three  types  of 
spectral  transition.  These  are 

(a)  bound-bound  transitions  which  give  rise  to  the  ordinary  line  spectra 
emitted  by  atoms  and  ions, 

(b)  free-bound  transitions  consisting  of  the  recombination  of  "free" 
plasma  electrons  into  bound  states  of  the  positive  ions  with  the  emission  of 
continuous  radiation,  and 

(c)  free-free  transitions  which  also  give  rise  to  continuous  radiation. 

This  process  can  be  described  classically  by  the  deceleration  of  electrons  in 
hyperbolic  orbits  in  the  fields  of  positive  ions,  thereby  transforming  kinetic 
energy  of  the  electron  into  radiation.  This  is  the  well  known  Bremsstrahlung 
radiation. 

In  view  of  the  extreme  conditions  of  temperature  and  pressure  in  the  plasma, 
line  spectra  will  be  broadened  and  shifted  by  perturbation  of  the  energy  levels 
by  particles  in  the  neighbourhood  of  the  radiating  atom.  Also,  when  there  is  an 
appreciable  degree  of  ionization,  the  electrical  fields  of  nearby  ions  will 
give  rise  to  Stark  broadening.  Trapped  magnetic  fields  in  the  plasma  may  also 
lead  to  Seeman  splitting.  Under  conditions  of  very  high  pressure  and 
temperature,  therefore,  most  of  the  line  radiation  will  be  smeared  out  to  such 
an  extent  that  it  is  masked  by  the  continuum.  Also,  when  the  degree  of 
ionisation  is  high,  the  occupation  number  for  a  particular  excited  state  is  low 
and  hence  the  line  it  emits  is  weak  and  merges  more  easily  with  the  continuous 
background.  Finally,  the  rare  gases  contain  a  large  number  of  closely  spaced 
lines  which  may  be  expected  to  overlap  under  these  conditions.  In  this  type  of 
plasma,  therefore,  the  contribution  of  bound-bound  transitions  to  the  total 
emitted  intensity  rail  be  small  and  is  neglected  in  the  intensity  calculations. 

In  the  case  of  optically  thick  plasmas,  in  local  thermodynamic  equilibrium, 
emitted  spectral  intensities  are  best  obtained  from  the  continuous  absorption 
coefficient  using  Kirchhoff’s  law  and  correcting  for  stimulated  emission. 
Allowance  must  also  be  made  for  absorption  of  radiation  by  the  plasma. 


/Lambert ' s 


-  7  - 

UNCLASSIFIED 


UNCLASSIFIED 


Lambert's  lav/  for  the  absorption  of  radiant  energy  is 

I  =  IQ  exp(-KA)  =  IQ  exp(-NoA)  .....  1j5 

where  I  is  the  intensity  of  the  incident  beam;  I  is  the  intensity  of  the  beam 
after  it  traverses  a  distance  Z  of  the  gas;  K  is  the  linear  absorption 
coefficient  of  the  gas;  a  is  the  cross-section  for  this  absorption;  and  N 
is  the  density  of  the  absorbing  particles.  N  and  a  must  be  individually 
considered  for  each  of  the  types  of  radiation  absorption  processes  and  if  there 
are  several  for  the  frequency  or  frequency  region  of  interest,  they  must  be 
summed. 

Assuming  that  the  concept  of  local  thermodynamic  equilibrium  is  valid,  or 
equivalently  that  Kirchhoff's  law  holds,  the  gaseous  emissitivity  of  the  pla3ma 
can  be  written  as  (lO) 

«■  =  1  -  exp(-KA)  . 14 

Following  Breene  and  Nardone  (i 1 )  a  factor  /3  is  introduced  into  the  exponent 
giving 

€  =  1  -  exp(-/?K&)  . 1j> 

In  deriving  the  above  emis3ivity  expressions,  one  must  consider  the 
radiation  arriving  at  the  surface  of  an  emitting  gaseous  layer  after  its 
augmentation  and  diminution  during  passage  through  the  layer  along  a  particular 
path.  At  a  given  point  on  the  boundary,  however,  radiation  is  arriving  along 
any  number  of  paths  and  an  integration  over  these  paths  is  necessary.  For 
paths  of  the  order  of  a  few  centimetres  and  smaller,  it  has  been  found  that  a 
value  of  /?  =  1.8  approximates  the  results  of  angular  integration. 

The  spectral  radiance  is  then  obtained  by  multiplication  with  the  Planck 
function 


Tv  =  evBvdv  =  Bv^  “  exp (-/?K&)]dv  . 1 6 

where  B(;dv  =  ^jpf—[exp(hv/kT)  -  l]_1dv 
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3*4  UnsUld-Kramers  Theory 


The  problem  is  now  resolved  into  one  of  calculating  the  continuous 
absorption  coefficient  Ky.  In  principle,  Ky  can  be  calculated  from  quantum 
mechanical  considerations,  but  to  the  present  time  this  has  only  been 
accomplished  for  one-electron  systems.  UnsiJld  (l2)  using  the  classical 
Kramers  (13)  theory,  has  developed  a  theoretical  expression  for  Bremsstrahlung 
and  recombination  radiation.  According  to  this  theory,  the  free-free 
absorption  index  of  radiation  of  frequency  v  in  the  field  of  a  singly  charged 
ion  is  given  by 


K  = 


16  7T 


2  _6 


N  2 
e 


3  ch(2jnn,P"  (kT)*  v 

3 


t  •  /  exp  (-a)  da 

2  v3  Jo 


. 17 


=  3.70  x  10 


N  ' 
,e  e 


.3 
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where  N  is  the  electron  density;  a  =  m  V2/kT,  V  is  the  electron  velocity  and 
©  © 

T  is  the  electron  temperature. 


If  local  thermodynamic  equilibrium  pertains,  then  Equation  1_8  can  be 
generalized  to  include  free-bound  absorption.  For  this  purpose,  the  lower 
integration  limit  should  not  be  zero  but  -a  =  -hv/kT,  corresponding  to  the 
limit  of  photoelectric  absorption  of  frequency  v.  The  total  absorption  index, 
taking  account  of  both  free-free  and  free-bound  absorption  of  a  single  charged 
ion,  will  have  the  form. 


K.  = 


V 


16  jrV  -  V  1  /°°  ,  V. 

— ■  "  - 2  •  8  *  . r  •  ~~~  /  exp  (-a;  da 

3  ./3  ch(2wm  )2  (kT)2  v 3  J- a 
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_  v  1 

=  3»70  x  10®  g  .  — r  *  —  exp(hv/kT) 

„3 


20 


In  the  case  of  an  optically  thick  plasma,  induced  emission  must  be 
considered.  The  absorption  coefficient  corrected  for  induced  emission,  can  be 
expressed  by 


K 


3.70  x  108  g 


A 

r  •  — [exp  (hv/kT)  -  1] 
T2  v3 


21 
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The  factor  g  has  been  introduced  to  correct  for  the  dissimilarity  between  the 
atom  and  hydrogen,  that  is,  partial  penetration  of  electrons  into  the  shell  of 
the  atom. 

The  applicability  of  this  equation  to  specific  atoms  can  be  evaluated  by 
considering  the  model  term  scheme  shown  in  Figure  2.  In  his  treatment,  Unstild 
arsumed  recombination  into  hydrogen-like  states  which  were  sufficiently  dense 
to  permit  the  sum  over  discrete  states  to  be  replaced  by  an  integral.  Thus 
Equation  22  gives  the  absorption  index  up  to  a  limiting  frequency  v^, 
corresponding  to  an  energy  hv  •  The  sequence  of  terms  is  considered  as  being 
sufficientl"'  dense  up  to  the  term  hi^  (measuring  from  the  ionization  limit),  to 
replace  term  summation  by  integration,  the  contribution  of  terms  below  hVg  must 
be  accounted  for  by  individual  summation.  If  no  terms  exist  between  hi4  and 
the  ground  state,  the  lower  boundary  of  the  integral  in  Equation  1 9  (for 
v  >  Vg)  will  remain  constant  and  equal  to  -a  =  -hv ^/kT.  In  ‘this  Rase,  the 
absorption  coefficient  will  be  given  by 

N  2  . 

Kv  =  3*70  x  108g  .-T-  .  —  [expChVg/kT)  -  l]  . 22 

-) 

and  the  absorption  index  will  decrease  according  to  -jjy.  A  comparison  of 
Equations  1_8  and  21_  gives  the  ratio  of  total  to  free-free  absorption 
coefficient  as 


g  [exp(hv/kT)  -  l] 


Assuming,  in  the  general  case,  that  g  =  1 ,  at  a  temperature  of  20,000°K, 
free-free  absorption  contributes  10  per  cent  of  the  total  at  3000  A  and 
43  per  cent  at  60C0  A.  At  9CC0  A,  the  contribution  of  Bremsstrahlung 
absorption  reaches  80  per  cent  of  the  total.  The  contribution  of  free-free 
absorption  is  thus  small  in  the  ultra-violet  region  but  increases  with  wave¬ 
length.  It  also  increases  in  importance  at  higher  temperatures. 
Bremsstrahlung  continuum  will  therefore  play  an  important  role  under  the 
conditions  of  interest  to  us. 

An  examination  of  the  energy  level  diagram  for  the  unperturbed  argon 
atom,  reveals  the  existence  of  a  continuous  sequence  of  terms  between  the 
ionisation  limit  and  the  3d6  term  at  v*  =  16,873  cm"1.  There  is  then  a 
window  of  the  order  of  3000  cm”1  in  the  term  sequence  until  the  next  series 
of  terms.  This  group  of  levels  starts  at  v*  =  19,818  cm"1  and  ends  at 
y*  =  24,439  cm"1  •  A  second  window  of  about  9CC0  cm"1  exists  between  the  2p1 0 
level  at  v?  =  ?4>439  cm"1  and  the  last  sequence  of  terras  starting  at 
V4  =  33,141  cm"1  and  ending  at  v|  =  35*397  cm"1.  There  are  no  terms  between 
v?  and  the  ground  level. 
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The  absorption  coefficient  should  therefore  be  determined  by  Equation 
up  to  a  limiting  frequency  (5  x  1014  sec”1),  corresponding  to  the  interval 
between  the  ionization  limit  and  the  3d6  term.  For  frequencies  between  vi  and 
i>2 }  the  absorption  index  will  vary  according  to  Equation  22.  To  account  for 
the  dense  sequence  of  terms  beginning  at  vt  and  ending  at  1/3,  the  absorption 
index  will  have  the  form  (leaving  out  the  induced  emission  term  for  clarity) 


-  H  5 

3.7  x  io8  .  -f 
v3  T2 


h  * 

V2 

kT 

e  + 


for  v  >  24,439  cm-1  to  d  =  33,141  cm”1  the  absorption  coefficient  will  be 
determined  by  the  expression 


T  N, 


K  =  3.70  X  108  4  .  -f-  [e 


kT-, 
+  e  ] 


for  the  frequency  region  from  V4  to  v5  the  index  will  have  the  form 


-  N  2  - 

K  =  3.70  x  108  4  .  -T-  [e 


h  * 

V2 

kT 

e  + 


kT-, 
+  e  J 


and  for  frequencies  v *  >  35,397  cm”1  the  appropriate  expression  will  be 


-  N_ 


=  3.70  x  10®  . 


3*5  Theoretical  Absorption  Coefficient  for  Argon 

The  calculated  absorption  index  for  argon  using  T  =  19,880°K, 

Ng  =  1*57  x  102°  cm”1  and  g  =  0.9,  is  shown  in  Figure  3»  It  is  compared  with 
the  absorption  coefficient,  calculated  assuming  Equation  22  is  valid,  down  to 
a  wavelength  of  2CC0  A.  This  assunption  was  made  by  Petschek  et  al.  (l4), 
who  noted  that  the  capture  of  electrons  with  energies  greater  than  about  */2  M 
to  the  first  excited  state,  will  be  fairly  infrequent.  The  capture  of  electrons 
to  the  ground  state  can  be  neglected.  The  first  excited  state  of  argon  is 


/about 


-  11  - 

UNCLASSIFIED 


UNCLASSIFIED 


about  4«2  eV  below  the  ionisation  limit  and  at  a  temperature  of  20,  COO0  K, 

3/2  kT  corresponds  to  about  2.6  eV.  This  leads  to  an  upper  frequency  cut-off 
which  corresponds  to  the  frequency  interval  Atf  given  by 


hAv  =  6.8  eV 


Thus  Equation  22  should  be  valid  for  a  frequency  interval  of  about  6.8  eV. 

In  high  density  ionised  plasmas,  however,  the  lowering  of  the  ionization 
limit  will  shift  the  cut-off  frequency  to  a  lower  value.  Since  accurate  values 
of  the  reduced  ionization  potential  are  not  available  at  present,  this  was  not 
taken  into  consideration.  This  could  lead  to  appreciable  error  when  large 
reductions  in  the  ionization  potential  exist. 

3.6  Theoretical  Absorption  Coefficient  for  Helium 


Using  Equation  22,  the  absorption  coefficient  for  helium  at  a  wavelength 
of  5CC0  A  was  calculated  for  incident  and  reflected  shocks.  Behind  the 
incident  sb'  ck  wave  (for  T  =  9465°K  and  =  2.63  x  1014  cm-3)  Ky  =  2.6  x  10“8 
cm”1,  and  behind  the  reflected  shock  (for  T  =  20,370°K  and  Ng  =  2.43  x  101 0 
cm”3),  a  0.29  cm”1.  These  represent  minimum  values  of  absorption 
coefficient  since  the  lowering  of  the  ionization  potential  was  not  taken  into 
account  in  calculating  the  temperature  and  electron  density.  The  effect  of 
the  fine  structure  of  the  helium  spectrum  on  the  absorption  coefficient  is  not 
evaluated  here. 

3.7  Calculated  Radiant  Emission 
(a)  Argon 

The  variation  of  emissivity  (calculated  from  Equation  15)  with  depth  of 
emitting  layer  at  a  wavelength  of  5CQP  A,  is  illustrated  in  Figure  4*  The 
optical  depth  (the  thickness  of  plasma  which  is  contributing  directly  to  the 
observed  spectral  intensity)  is  about  0.C6  mm.  When  Equation  16  is  valid, 
it  is  apparent  that,  if  the  thickness  of  the  emitting  layer  is  greater  than  the 
optical  depth  at  a  given  wavelength,  then  the  plasma  should  radiate  as  a  black 
body  at  the  equilibrium  temperature.  It  also  follows  that,  if  energy 
relaxation  times  are  negligible  behind  the  shock  front,  the  rise  time  of  the 
radiation  at  5CC0  A  should  be  about  10”8  sec  and  the  continuum  intensity 
profile  should  resemble  that  of  the  emissivity  curve  shown  in  Figure  4*  The 
ionization  relaxation  times  behind  the  shock  front  were  calculated  from  the 
expression  ' 


P  =  0.156  exp(87,000/T) 

/where 
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where  P  is  the  pressure  in  torr  and  t,  is  the  ionization  relaxation  time  in 
microseconds.  Using  this  equation,  was  found  to  be  1.6  x  10“8  sec.  This 
expression  was  obtained  by  Petschek  and  Byron  (15)  using  high  purity  argon 
and  their  results  indicate  a  marked  decrease  in  relaxation  times  when  less 
pure  gas  was  employed.  The  calculated  relaxation  time,  therefore,  represents 
a  maximum  value  and  the  actual  time  is  probably  considerably  less.  Electrons 
ejected  into  the  argon  from  the  detonating  explosive  will  also  reduce  the 
relaxation  time.  Energy  relaxation  times  are  therefore  probably  negligible 
compared  with  the  radiation  rise  times. 

(b)  Helium 

Since  the  helium  behind  the  incident  shock  wave  has  a  very  low  absorption 
coefficient,  its  emissivity  will  be  low  and  very  little  radiation  will  be 
observed  for  emitting  layers  of  a  few  centimetres  and  less.  The  variation  of 
emissivity  with  thickness  of  helium  plasma  at  50C0  A  behind  the  reflected  shock 
wave,  is  shown  in’ Figure  5-  The  optical  depth  is  about  10  cm,  but  this  must  be 
regarded  as  a  maximum  value. 


4.  SPECTRAL  MEASUREMENTS 


4.1  Spectroscopic  Techniques 

The  absolute  intensity  and  spectral  distribution  of  the  energy  radiated  by 
shock-heated  gases  were  measured  using  a  Hilger  medium  quartz  spectrograph 
equipped  with  photoelectric  recording  facilities.  The  intensity  of  the 
radiation  from  the  shock-heated  plasma  was  compared  at  six  selected  narrow 
wavebands  in  the  visible  spectral  region  with  that  of  a  tungsten  ribbon- 
filament  lamp,  previously  calibrated  by  the  National  Physical  Laboratory  for 
brightness  temperature  against  current.  Identical  optics  were  used  for  each  of 
the  two  radiators  in  each  experiment,  so  that,  in  calibrating  the  photo¬ 
multipliers,  the  spectrograph  and  the  associated  optics  were  also  calibrated. 
Reflective  optics  were  used  to  focus  the  radiation  on  to  the  entrance  slit  of 
the  spectrograph.  The  spectral  intensity  of  the  lamp  at  the  desired  wavelength 
and  temperature  was  calculated  employing  data  on  the  emissivity  of  tungsten 
given  by  De  Vos  (16).  The  output  of  each  photomultiplier  was  displayed  as  a 
function  of  time  on  an  oscilloscope  and  recorded  photographically,  the  overall 
rise- time  of  the  equipment  being  15  nanosec.  To  achieve  a  wide  range  of  light 
sensitivities,  the  photomultiplier  anode  loads  were  varied.  Since  the  photo¬ 
multipliers  were  calibrated  at  a  continuous  low  output  current  level  and  then 
used  to  measure  high  current  pulses,  care  was  taken  to  operate  the  tubes  only 
on  the  linear  portions  of  their  response  curves.  There  are  several  possible 
sources  of  error  in  the  measurements,  such  as  error  in  ribbon  lamp  brightness 
•  tempei'ature  (±  7  aegC),  short-term  drift  of  the  photomultiplier  sensitivities, 
uncertainties  in  the  standard  lamp  current  measurement  and  in  the  positioning 
of  the  lamp,  so  that  a  precise  determination  of  error  is  difficult.  However, 
the  maximum  error  in  intensity  measurement  is  probably  less  than  15  per  cent. 

/4.2  . 
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4*2  Continuum  Intensity  ileasurements 


The  radiative  properties  of  argon  heated  by  shock  waves  generated  by 
various  detonating  high  explosives  were  investigated.  The  hot  plasma  was 
viewed  normal  to  the  direction  of  the  shock  frorit.  Radiation  reaching  the 
spectrometer  originated  from  only  a  small  area  about  2  mm  x  0.2  mm  at  the 
centre  of  the  shock  front.  Interference  from  edge  effects  arising  from 
interaction  of  the  shock  wave  with  the  walls  of  the  container  were  thus 
minimised.  The  effects  of  the  high  explosives  used  to  drive  the  shock  on  the 
continuum  intensity  was  examined  in  the  apparatus  shown  in  Figure  6.  The 
results  obtained  using  argon  at  1  atm  pressure  are  shewn  in  Figure  8.  In  the 
visible  spectrum  the  radiation  is  similar  to  that  from  a  black  body  at  about 
21,CC0°K.  Little  significant  difference  in  intensity  or  spectral  distribution 
was  found  using  luiX,  RT  (a  cast  explosive  containing  about  60  per  cent  RDX  and 
40  per  cent  TNT) ,  tetryl  or  a  tetryl/TNT  mixture  although  HLiX  gave  on  the 
average  more  light  than  RT.  Several  experiments  were  generally  performed  for 
each  explosive  composition.  There  was  a  variation  of  up  to  20  per  cent  in  the 
intensity  using  the  same  explosive.  Part  of  this  is  due  to  experimental  error 
but  variations  in  the  explosive  charge  such  as  composition,  uniformity  and 
physical  condition  probably  also  contribute. 

The  radiation  from  different  shocked  noble  gases  wa3  compared  using  the 
experimental  arrangement  illustrated  in  Figure  7*  Spectroscopic  grade  gases 
were  sealed  into  the  glass  tubes  at  a  pressure  of  440  torr,  using  high  vacuum 
techniques,  and  then  shocked  by  charges  of  RT.  The  impurity  level  was  less 
than  10  ppm  for  A,  He  and  He  while  Kr  and  Xe  were  99*5  per  cent  pure,  the 
balance  being  Xe  and  Kr  respectively.  In  the  case  of  argon,  krypton  and  xenon, 
intense  radiation  was  emitted  from  behind  the  incident  shock  waves,  with  about 
a  50  per  cent  increase  in  intensity  on  reflection  at  the  end  window.  However, 
the  rise  times  of  the  radiation  varied  considerably,  decreasing  in  the  order 
argon  (7*5  psec),  krypton  (l  psec)  and  xenon  (0.1  psec).  The  duration  of  the 
reflected  shock  also  decreased  in  the  order  of  increasing  molecular  weight 
from  helium  (2  psec)  to  xenon  (0.1  psec).  Typical  records  are  illustrated  in 
Figures  9  and  10.  No  light  emission  was  observed  from  behind  the  incident 
3hock  waves  in  neon  and  helium,  but  on  reflection  at  the  end  window  of  the 
tube  a  pulse  of  radiation  comparable  in  intensity  to  that  obtained  from  argon 
was  observed.  The  observed  radiancies,  which  are  compared  in  Figure  11, 
differed  by  less  than  a  factor  of  two  (two  experiments  were  performed  for  each 
gas). 


Experiments  were  performed  with  shock  waves  generated  by  the  explosion  of 
RT  in  cylinder  quality  helium,  using  an  apparatus  similar  to  that  shown  in 
Figure  6.  The  pathlength  of  helium  was  25  cm.  Very  little  radiation  was 
emitted  from  behina  the  incident  shock,  but  on  reflection  a  sharp  pulse  of 
radiation  lasting  for  about  2.5  p sec  was  obtained.  A  typical  record  at  two 
of  the  wavelengths  is  shown  in  Figure  12.  It  is  apparent  from  the  record 
that  the  radiation  is  cut  off  by  the  contact  surface  before  the  maximum 
intensity  is  attained. 
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4.3  Streak  Photographs  and  Shock  Velocities 

A  Beckman  and  1/hitley  l.'cdel  339  streak  camera  was  used  to  measure  the 
velocity  of  the  incident  shock  naves  and  of  the  contact  surface  between  the 
explosive  and  the  gas.  The  depth  of  the  editing  layer  behind  the  incident 
and  reflected  shocks  was  also  measured.  The  streak  photograph  of  73 
detonating  into  atmospheric  pressure  argon.  Figure  13A,  shows  the  distances 
of  the  shock  front  and  contact  surfaces  as  a  function  of  time.  The  streak 
camera  slit  was  oriented  parallel  to  the  tube  axis.  Figure  14  sets  out  the 
interpretation  of  the  streak  photographs  obtained.  The  incident  shock  wave 
propagates  down  the  column  of  argon  at  almost  constant  velor'ty  until  it  meets 
the  end  window  of  the  tube  when  it  is  reflected  back  into  *  i  hot  gas.  Ko 
increase  in  intensity  is  detectable  when  the  incident  shock  is  reflected.  The 
thickness  of  the  emitting  gas  layer  is  determined  by  the  velocity  of  the 
contact  surface  which  rapidly  cools  the  hot  gas,  quenching  the  radiation.  The 
maximum  thickness  of  hot  plasma  'at  the  end  o f  the  tube  was  found  to  be 
1.5  cm,  whilst  the  calculated  value  for  this  length  of  tube  was  1.2  cm  (using 
the  particle  velocity  computed  in  Section  3)»  If  this  difference  is 
significant,  it  may  mean  that  the  hot  gas  is  not  quenched  immediately  by  the 
explosive  products.  The  thickness  of  the  reflected  shock-heated  gas  will  of 
course  depend  on  the  relative  velocities  of  the  reflected  shock  wave  and  the 
contact  surface. 

’with  helium  (Figure  133)  no  radiation  was  observed  frra  the  incident 
shock  front  but  the  gas  behind  the  reflected  shock  was  highly  luminous.  A 
fairly  sharp  quenching  was  observed  when  the  contact  surface  reached  the 
radiating  plasma  but  the  gas  continued  to  radiate  at  a  much  lower  level  for 
some  time.  Ko  satisfactory  explanation  of  this  phenomenon  has  yet  been 
formulated. 


5.  DISCUSSION 


It  was  shown  in  Section  3  that  if  the  mean  free  path  for  radiation  of 
a  given  frequency  was  less  than  the  depth  of  the  emitting  region  the  emitted 
intensity  will  be  that  of  a  black  body.  From  the  very  rapid  rise  times  of  the 
radiation  in  argon  shocked  by  BT  (<0.1  psec)  the  optical  depth  must  be  scalier 
than  0.5  cm  and  is  probably  considerably  less.  Furthermore,  since  the  plasma 
achieves  a  depth  of  about  15  mm  after  traversing  15  cm  of  gas  it  should 
radiate  according  to  the  black  body  law.  The  observed  spectral  distribution 
(Figure  8)  closely  approximates  to  that  of  a  black  body  at  2l,0C0°X. 

In  view  of  the  high  electron  concentration  in  the  argon  plasma,  absorption 
of  radiation  will  be  very  high.  This  results  in  a  much  greater  cooling  time 
for  the  plasma  than  if  it  were  optically  thin.  Assuming  black  body  radiation, 
the  natural  radiative  cooling  time  of  the  argon  plasma  is  estimated  to  be 
greater  than  5C0  psec.  Since  the  plasma  is  only  allowed  to  radiate  for  a 
shc.r);  time  before  i z  is  quenched  by  the  contact  surface,  which  follows  the 
shock  front  very  closely,  this  represents  a  major  energy  loss  process  ana 
source  of  inefficiency. 
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To  obtain  a  high  radiance  a  high  temperature  is  necessary,  but  increasing 
the  shock  velocity  in  argon  will  not  have  much  effect  on  radiance  since  this 
will  mainly  result  in  increased  ionization  and  a  smaller  optical  depth.  In 
effect,  the  argon  temperature  is  buffered  by  thermal  ionization.  However,  if 
the  shock  velocity  is  sufficient  to  ionize  most  of  the  argon  atoms,  then  an 
abrupt  rise  in  temperature  should  occur  before  second  ionisation  becomes 
appreciable.  It  would  appear  preferable  at  first  sight,  to  employ  a  gas  of 
high  ionization  potential  such  as  helium  but  unfortunately  very  high  shock 
velocities  are  necessary  with  helium  in  order  to  produce  a  Mach  number 
comparable  to  that  in  argon.  Neon  might  prove  to  be  a  suitable  compromise  gas. 

An  alternative  approach  would  be  to  choose  a  material  of  very  low  first 
ionization  potential  such  as  caesium  which  is  completely  ionized  at  comparatively 
low  temperatures.  Since  the  second  ionization  potential  is  high  (23.4  eV)  it 
should  be  possible  to  produce  a  high  temperature,  optically  dense  plasma  with 
moderate  shock  velocities.  Unfortunately,  there  are  severe  practical  problems 
associated  with  the  use  of  caesium  vapour. 

High  temperature  black  body  sources  emit  most  ol  their  radiant  energy  in 
the  extreme  ultra-violet  region  but  this  is  of  no  use  for  pumping  the  majority 
of  laser  materials.  The  conversion  of  this  energy  into  visible  radiation 
would  result  in  a  major  improvement  in  the  efficiency  of  the  source. 

Experiments  with  fluorescent  convertors  have  so  far  proved  inconclusive.  A 
more  promising  approach  would  be  to  utilize  the  absorbing  characteristics  of 
the  gas  ahead  of  the  shock  wave  to  tailor  the  spectral  distribution  to  that  of 
the  pumping  band  of  a  given  laser,  perhaps  by  the  addition  of  small  amounts  of 
additives.  In  this  way  it  might  be  possible  to  achieve  a  non-equilibrium 
spectral  distribution. 
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APPENDIX 

Symbols 

temperature,  °K 
density,  g  cm"3 
velocity,  cm  sec"1 
atomic  mass,  g 

electron  mass,  9*1072  x  10"28  g 

degree  of  ionization 

ionization  potential,  eV 

gas  constant  per  unit  mass,  ergs  deg"1  g”1 

Boltzmann's  constant,  1*3803  x  10"16  erg  deg"1 

charge  of  electron,  4*803  x  10”10  eou 

electronic 

Planck's  constant,  6.6237  x  10"27  erg  sec 
partition  function 
pressure,  dyn  cm"2 
statistical  weight 

velocity  of  light,  2*9979  x  1010  cm  sec"1 
electron  density,  cm"3 
intensity,  erg  cm"2 
absorption  coefficient,  cm"1 
absorption  cross-section,  era”2 
emissivity 

/Symbols  (Cont'd.) 
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APPENDIX 

Symbols  (Cont'd. ) 

v  frequency,  sec**1 

-&  thickness  of  gas  layer,  cm 

th  , 

n  n  term 

N  density  of  absorbing  gas  particles,  cm"3 

j3  correction  term  for  emissivity 

g  correction  factor  for  absorption  coefficient 

3^  Planck  black  body  function,  erg  cm"2 
r,  ionization  relaxation  time,  10"6  sec 

1  gas  behind  shock  wave 

A  argon 

Kr  krypton 

Xe  xenon 

He  helium 

Ne  neon 

Subscripts 

o  undisturbed  gas 

1  gas  behind  shock  wave 
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FIG.  G.  SCHEMATIC  DIAGRAM  OF  SHOT  ASSEMBLY 
USING  .ATMOSPHERIC  PRESSURE  GASES. 
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USING  SPECTROSCOPIC  GASSES  AT  440  TORR. 


UNCLASSIFIED 


(_porsr}5  £_uj3  non  7jjsN3±NI  "WaiDBdS 


UNCLASSIFIED 


FIG.  8.  SPECTRAL  INTENSITY  FOR  ARGON  SHOCKED  8V  RDX/TNT 
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The  visible  radiation  from  a  thermally  Ionized  high  density  plasma 
I  produced  by  explosively  generated  shock  waves  has  been  Investigated  os  0 
!  possible  loser  pumping  source.  It  Is  considered  that  the  major  problem  Is  ■ 
!  the  efficiency  of  conversion  from  chemical  energy  to  radiation  In  the 
i  absorption  bond  of  the  loser.  The  study  has  therefore  been  concentrated  on 
|  the  mechanism  of  the  emissive  process  with  a  view  to  increasing  the  tudlution 
j  emitted  In  a  given  spectral  region.  The  real  cos  thermodynamic  variables 
i  behind  the  shock  waves  have  been  evaluated  for  argon  and  helium  end  used  to 
j  calculate  linear  absorption  coefficients  based  on  UnsBld-Kremers  theory, 

:  Using  Klrchhoff's  low,  the  radiant  Intensities  at  the  surfaces  of  optically 
J  thick  argon  and  helium  plasmas  at  20,CC0°K  have  been  calculated  and  compared  j 
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No.  12/R/67 

R.T.  Bailey  and  J.T.A.  Burton  November  1967 


The  visible  radiation  rrem  a  thermally  Ionized  high  density  plasma 
produced  by  explosively  generated  shock  waves  has  been  Investigated  as  a 
possible  laser  pumping  source.  It  Is  considered  that  the  major  problem  is 
the  efficiency  of  conversion  from  chemical  energy  to  radiation  in  the 
absorption  band  of  the  laser.  The  study  has  therefore  been  concentrated  on 
the  mechanism  of  the  emissive  process  with  a  view  to  increasing  the  radiation 
emitted  in  a  given  spectral  region.  The  real  gas  thermodynamic  variables 
behind  the  shock  waves  have  been  evaluated  for  argon  and  helium  and  used  to 
calculate  linear  absorption  coefficients  based  on  U’nsOld-Krar.ers  theory. 

Using  Klrchhoff's  law,  the  radiant  Intensities  at  the  surfaces  of  optically 
thick  argon  and  helium  plasmas  at  20,C0C°K  have  been  calculated  and  compared 
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